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1.  Introduction 


The  37-mm  chamber  (figure  1)  is  a  proven  test  bed  for  evaluating  propellant  effects  in  realistic 
ballistic  environments.  The  U.S.  Army  Research  Laboratory  (ARL)  has  been  tasked  to  evaluate 
the  performance  of  ceramic  nozzles  in  these  ballistic  environments.  The  area  of  the  gun  that 
appears  to  wear  rapidly  is  near  the  origin  of  rifling  where  obturation  is  poor.  Ceramics  offer 
particularly  good  heat  resistance  and  compressive  strength  and  are  therefore  a  logical  choice  to 
reduce  wear  in  these  areas  and  prolong  life  of  the  barrel.  This  is  an  area  of  concern  as  the 
high-performance  propellants  now  in  use  for  tactical  projectiles  can  result  in  a  condemned  tube 
in  as  few  as  175  rounds.  Ceramic  inserts  are  therefore  certainly  worthy  of  examination. 


Figure  1.  ARL  37-mm  erosion  fixture. 

One  of  the  requirements  for  a  successful  insert  is  its  ability  to  withstand  not  only  high 
temperatures  but  also  high  pressures  simultaneously.  A  series  of  nozzles  with  different  ceramic 
inserts  were  fabricated  and  placed  in  the  37-mm  erosion  fixture.  The  results  observed  in  many  of 
the  samples  suggested  that  examination  of  the  stress  states  during  firing  was  warranted  and 
would  help  indicate  sources  of  the  failure.  These  stresses  are  caused  by  both  thermal  and 
pressure  loads.  A  finite-element  model  of  the  37-mm  chamber  was  created  to  predict  the 
thermal  evolution  and  corroborate  the  failure  modes  observed. 

In  this  analysis,  thermal  conductive  theory  and  modeling  methods  are  presented  first.  The 
finite-element  thermal  analysis  is  then  performed  for  the  full-scale  gun  chamber-nozzle  model. 

A  space-  and  time-dependent  convective  boundary  condition  is  applied  to  the  chamber  and 
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nozzle.  The  peak  temperatures  in  the  ceramic  nozzles  and  their  temperature  gradients  are 
determined  and  discussed. 

Figure  2  illustrates  a  quarter  of  a  three-dimensional  (3-D)  finite-element  model  for  the  37-mm 
gun  test  fixture.  The  chamber  is  fixed  on  a  solid  base,  and  the  nozzle  is  placed  at  the  front  of  the 
chamber  backed  up  with  a  rupture  diaphragm.  A  steel  retainer  is  used  to  tighten  the  system,  as 
shown  in  the  figure  2.  As  the  propellant  reaches  burnout ,  the  high-pressure  ruptures  the 
diaphragm. 


Three  candidate  ceramic  materials  are  selected  in  this  investigation.  SN47  silicon  nitride 
(Si3N4),  produced  by  Ceradyne,  Inc.,  is  a  sintered  reaction  bonded  material  that  contains  a  small 
amount  of  yttria  (Y203)  and  alumina  (A1203)  as  sintering  aids.  The  fabrication  process  results 
in  a  microstructure  containing  interlocking  needle-like  grains.  STK4  is  a  SiAlON  material 
produced  by  Kennametal  using  sintering  and  a  post-sintering  hot  isostatic  pressing  step  to 
produce  a  fully  dense  material.  The  microstructure  contains  elongated  grains  and  consists  of 
60%  alpha  phase  and  40%  of  the  beta  phase.  ZR02  is  sintered  zirconia  (Zr02)  produced  by 
CoorsTek  using  ceria  (Ce02)  to  retain  the  meta-stability  and  the  tetragonal  phase  at  room 
temperature.  This  results  in  a  transformation  toughened  material  with  high  strength  and 
toughness  at  room  temperature. 

In  the  simulation,  the  steel  nozzle  is  modeled  as  the  baseline.  Then,  the  previously  mentioned 
three  ceramic  nozzles  are  simulated  and  compared  with  the  baseline. 


2.  Theory  Background  and  the  Chamber  Gas  Temperature  Condition 


Application  of  the  first  law  of  thermodynamics  to  the  control  volume  yields  the  general 
conduction  equation  for  homogeneous  isotropic  solid  materials  (7): 
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(1) 


d  (1dT{x,y,z,t)  Q  dT(x,y,z,t)  Q  dT(x,y,z,t )  dT(x,y,z,t) 

dxKK  )  +  ZW  )  +  ^  >-PLp  yrt  • 

The  temperature  T  is  defined  in  the  domain  (x,  y,  z)  as  a  function  of  time,  t.  Here,  p  is  the 
density,  and  Cp  is  the  specific  heat  per  unit  mass.  Isotropic  thermal  conductivity  is  represented 
by  K.  In  the  current  case,  both  specific  heat  Cp  and  thermal  conductivity  K  are  temperature 
dependent.  Since  the  solution  of  the  differential  equation  1  involves  a  number  of  integration 
constants,  the  completion  of  the  formulation  requires  an  equal  number  of  appropriate  boundary 
conditions  in  space  and  time  to  determine  these  constants. 

2.1  Initial  and  Boundary  Conditions 

For  an  unsteady  problem,  the  temperature  of  the  continuum  under  consideration  must  be  known 
at  some  instant  of  time.  In  many  cases,  this  instant  is  most  conveniently  taken  to  be  the 
beginning  of  the  problem.  If  the  initial  condition  is  given  by  To,  the  solution  of  this  problem 
returns  T (r,z,t)  at  all  points  of  the  continuum  to  be  specified  : 

T(r,z,t=0)  =  T0 .  (2) 

To  is  taken  to  be  the  ambient  temperature  in  our  simulation. 

The  most  commonly  encountered  boundary  conditions  are  prescribed  temperature,  insulated 
wall,  and  convective  heat  transfer.  These  three  boundary  conditions  can  be  unified  as  one 
generalized  temperature  boundary  condition  mathematically  as: 

dT 

a  — -  =  h  [bT^  -  cT  (?))  (3) 

on 

where  o/an  denotes  the  differentiation  along  the  normal,  Ts  is  the  surface  temperature,  and  T(t)  is 
the  temperature  at  a  distance  far  from  the  boundaries.  The  term,  h,  defines  the  effective  heat 
transfer  coefficient  quantifying  the  convective  heat  flux  between  the  surface  and  the 
environment,  and  the  constants  a,  b,  c  are  either  1  or  0  depending  on  the  type  of  boundary  heat 
transfer.  If  the  surface  temperature  of  the  boundaries  is  specified  to  be  a  function  of  space  or 
time,  the  coefficient  a  is  zero  and  b  and  c  are  equal  to  one.  When  the  heat  transfer  across  the 
boundaries  of  a  continuum  cannot  be  prescribed,  the  constants  ( a ,  b,  and  c)  are  all  equal  to  one. 
The  heat  flux  across  the  surface  is  then  proportional  to  the  temperature  difference  between  the 
boundary  and  the  driving  temperature.  In  the  current  case,  the  outside  chamber  wall  is  ambient 
temperature.  The  inside  wall  heat  transfer  boundary  condition  is  convective  and  is  discussed  in 
the  next  section. 

2.2  Gas  Temperature  and  the  Chamber  Boundary  Condition 

In  the  past,  extensive  interior  ballistic  (IB)  tests  and  calculations  were  conducted  to  determine 
the  chamber  temperature  and  internal  pressure  for  the  several  propellants  (2).  Inside  the 
chamber,  the  gas  temperature  was  found  to  be  space  dependent.  Especially  in  the  nozzle  area, 
the  gas  temperature  has  a  noticeable  change  due  to  the  significant  geometric  variation.  IB 
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calculations  coupled  with  a  finite  difference  heat  conduction  code  XBR2D  v54  (REF  Conroy) 
produced  the  gas-temperature  profiles  as  well  as  heat-transfer  coefficients  as  a  function  of  time 
and  axial  location.  Figure  3  shows  the  ex99  propellant  gas-temperature  profiles  in  the  chamber. 
It  is  noted  that  a  200  K  difference  occurred  for  the  peak  temperature  in  axial  distance  from  226.1 
to  238.8  mm,  which  is  the  exact  nozzle  position  (see  figure  4). 


Figure  3.  37-mm  erosion  fixture  gas-temperature  profiles. 


The  convective  coefficient  h  detennines  the  driving  gas-temperature  difference.  In  a  gun  system, 
the  convective  coefficient  is  no  longer  constant  and  behaves  dynamically.  In  the  current  case, 
the  convective  coefficient  is  time  and  space  dependent.  According  to  the  experimental  data,  the 
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coefficient  profile  was  obtained  using  the  ARL  Gun  Tube  Heat  Transfer  Code  (3)  (2003 
version),  as  shown  in  the  figure  5.  The  coefficient  is  zero  before  the  diaphragm  ruptures  at 
12.5  ps.  The  corresponding  chamber  inside  surface  is  insulated  from  the  gas  temperature.  The 
peak  values  occur  at  12.5  ps,  after  which  they  sharply  decrease.  In  the  nozzle  area,  the 
coefficients  vary  as  much  as  38%  in  the  axial  direction  due  to  the  flow  conditions. 


Figure  5.  Time-  and  space-dependent,  heat-convective  coefficient  profiles. 


3.  Finite-Element  Approach  (FEA) 


When  the  boundaries  are  irregular,  or  when  the  problems  involve  multiple  domains,  there  is  no 
exact  analytical  solution  available  for  the  heat-conduction  problem  described  by  equation  1 . 
Approximate  integral  and  variational  techniques  are  not  flexible  enough  to  handle  the 
multidomain  problem.  In  the  past,  the  finite-difference  method  was  widely  used  to  solve  such 
heat-conducting  problems.  However,  for  3-D  multidomain  problems,  the  finite-difference 
method  is  extremely  tedious  and  time  consuming  ( 4 ).  In  contrast,  the  finite-element  method 
offers  promising  advantages.  The  finite-element  code  ABAQUS,  developed  by  Hibbitt,  Karlsson 
&  Sorensen,  Inc.,  is  used  to  solve  the  heat-conducting  problem  formulated  above  (5). 

Using  ABAQUS,  there  is  no  limitation  to  the  complexity  of  the  geometry  or  the  number  of 
domains.  The  structure  can  be  meshed  using  any  suitable  pre-processor  code  or  directly  using 
ABAQUS  meshing  input  cards  according  to  the  user’s  menu,  as  was  done  for  this  analysis  ( 6 ). 
The  space-  and  time-dependent,  gas-temperature  input  and  convective  coefficient  h  for  the 
fixture  can  be  coded  by  using  user  subroutines  interfaced  with  the  ABAQUS  main  code.  The 
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formulations  for  the  temperature  and  convective  coefficient  shown  in  figures  2  and  3  are 
externally  programmed  and  linked  into  existing  ABAQUS  code. 

In  the  model,  shown  in  figure  5,  the  interfaces  along  some  of  the  different  components  are  fixed 
and  the  adjoining  elements  share  common  nodes  at  the  interfaces.  The  total  chamber  length  is 
238.8  mm  and  the  nozzle  length  is  12.7  mm.  The  chamber  wall  thickness  is  25.4  mm.  Three- 
dimensional  eight-node  linear  brick  elements  were  used  in  the  model.  Extra  care  was  taken  for 
meshing  each  part  of  the  model,  especially  along  their  common  edges  and  nozzle  area.  A  finer 
mesh  is  used  for  the  nozzle  and  its  vicinity.  Only  one  quarter  of  the  model  is  used  in  the  finite- 
element  model  simulation  with  the  symmetric  planes  and  the  two  ends  having  insulated 
boundary  conditions  assigned.  The  ABAQUS  CAE  software  is  used  for  post  processing  which 
includes  contour  plots  and  temperature  profile  curves. 

3.1  Input  Material  Properties 

In  the  simulation,  the  temperature  varies  from  room  temperature  to  a  peak  of  1200  °C.  Over 
such  a  large  temperature  range,  the  material  thermal  conductivities  and  heat  capacities  are  no 
longer  constants.  The  thermal  conductivities  of  chamber  and  ceramic  nozzles  are  plotted  in 
figure  6.  It  can  be  seen  that  the  steel  thennal  conductivity  changes  much  more  dramatically  with 
the  temperature  than  that  of  the  ceramic.  The  thermal  conductivities  of  ceramics  ZR02  and 
STK4  vary  little  as  the  temperature  is  elevated  and  therefore  they  are  treated  as  constants  in  the 
simulation. 


Figure  6.  Thermal  conductivities  of  steel  4340  and  ceramics  (7,  8). 
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A  similar  situation  is  found  for  the  heat  capacity  of  these  materials  and  is  shown  in  figure  7.  All 
three  ceramics  have  thermally  varying  heat  capacities  from  room  temperature  to  1000  °C. 

The  material  density  and  coefficient  of  thermal  expansion  (CTE)  for  steel  4340  and  ceramics  are 
listed  in  table  1 .  It  is  noted  that  SN47  and  STK4  are  very  close  in  their  densities  and  CTE. 


Temperature  (C) 


Figure  7.  Heat  capacities  of  steel  4340  and  ceramics  (7). 


Table  1 .  Density  and  coefficient  of  thermal  expansion  of 
steel  4340  and  ceramics  (7). 


Materials 

Steel  4340 

SN47 

STK4 

ZR02 

Density  (kg/m3) 

7800 

3200 

3400 

6100 

CTE  (  10_6/°C) 

8.4 

3.2 

3.3 

11.8 

3.2  Model  Convergence  Study 

The  finite-element  analysis  simulation  time  step  was  set  at  a  0.1 -jus  increment  over  a  total  time 
duration  of  3 1  ps.  The  number  of  elements  used  was  based  on  the  model  convergence  study. 
Figure  8  shows  the  model  convergence  result  with  number  of  elements  vs.  the  normalized  nozzle 
temperature.  There  were  a  total  of  10  runs  for  the  different  model  sizes  from  4000  to  14,000 
elements.  A  normalized  nozzle-edge  temperature  is  defined  as  the  ratio  of  the  nozzle-edge 
temperature  of  the  model  with  various  elements  versus  the  temperature  at  the  same  location  with 
14,000  elements.  Using  9000  elements,  the  model  results  converge.  Further,  if  only  the  nozzle 
area  was  of  interest,  even  fewer  elements  could  be  used  by  reducing  number  of  elements  in  the 
chamber  and  other  parts. 
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Figure  8.  The  variation  of  the  nozzle  temperature  with  the  number  of  elements. 


4.  Heat  Conductive  Analysis  of  the  37-mm  Gun  Chamber 


The  simulation  of  the  steel_4340  gun  chamber  with  a  steel_4340  nozzle  is  used  as  a  baseline  for 
comparison  to  nozzle  ceramic  materials  SN47,  STK4,  and  ZR02.  For  the  convenience  of 
analysis,  areas  of  interest  in  the  nozzles  are  numbered  as  shown  in  figure  9.  The  nozzle  entrance, 
throat,  and  top  areas  are  denoted  as  number  1,  2,  and  3,  respectively.  The  gun  chamber  is 
denoted  as  number  4. 

4.1  Temperature  Contour  for  the  Steel  Chamber  With  Steel  Nozzle 

Figure  9a  shows  the  temperature  distribution  in  the  steel_4340  nozzle  and  its  vicinity  at  the  IB 
time  17  ms,  the  time  of  peak  temperature  in  the  chamber.  It  is  seen  that  the  nozzle  entrance 
(area  denoted  1)  has  the  highest  temperature  of  1650  °C.  In  the  nozzle  top  area  (area  3),  the 
temperature  is  about  900  °C,  while  in  the  throat  (area  2)  the  temperature  reaches  820  °C. 

In  the  chamber  region  (area  4),  the  temperature  is  approximate  550  °C.  The  temperature  gradient 
is  captured  in  the  simulation  along  the  nozzle  axial  direction  primarily  because  of  the  space- 
dependent  gas  temperature  and  convective  coefficient. 
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Figure  9.  Temperature  contours  at  the  highest  levels  for  the  steel  and  ceramic  nozzles. 

4.2  Temperature  Contours  for  the  Steel  Chamber  With  Ceramic  Nozzles 

The  peak  temperature  distributions  in  the  ceramic  nozzles  with  the  steel_4340  gun  chamber  are 
shown  in  figure  9b  for  SN47,  figure  9c  for  STK4,  and  figure  9d  for  ZR04.  The  thermal  response 
not  only  depends  on  the  material  conductivity,  but  also  on  the  material  density  and  heat  capacity 
(see  equation  1).  For  the  SN47  ceramic  nozzle  (figure  9b),  the  temperature  in  the  nozzle 
edge  can  be  as  high  as  1850  °C.  In  the  nozzle  throat  and  top  areas,  the  temperatures  are  about 
1200  and  1150  °C,  respectively.  The  SN47  nozzle  is  about  250  °C  higher  than  that  computed  for 
the  steel_4340  nozzle.  In  the  chamber  area  (area  4),  the  temperatures  are  all  the  same,  regardless 
of  the  nozzle  material.  The  temperature  contour  for  STK4  nozzle  (figure  9c)  is  quite  close  to  that 
of  SN47.  The  STK4  nozzle  is  about  30  °C  lower  than  that  of  SN47  at  the  nozzle  entrance 
location.  By  reviewing  the  thermal  properties  of  these  two  ceramics,  it  is  noted  that  they  are 
quite  similar  in  density,  conductivity,  and  heat  capacity  (figures  6  and  7  and  table  1).  It  therefore 
makes  sense  that  they  have  quite  similar  thermal  response.  For  the  ZR02  ceramic  nozzle,  the 
temperature  is  about  280  °C  higher  than  that  of  the  steel  nozzle  at  the  nozzle  edge,  throat,  and  top 
areas. 

4.3  Temperature  Profiles 

The  previous  discussion  is  based  on  the  peak  temperature  response.  It  is  beneficial  to  study  the 
temperature  development  throughout  the  ballistic  cycle.  Figure  10  shows  the  temperature 
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profiles  as  a  function  of  time  for  the  steel_4340  baseline  and  three  ceramic  nozzle  materials  in 
the  nozzle  edge  and  chamber  areas.  Corresponding  to  the  thermal  convective  coefficient  curve 
in  figure  5,  the  temperatures  remain  ambient  until  the  diaphragm  ruptures  at  12.5  ps.  Within 
4  ps,  the  temperatures  sharply  increase  to  their  peak  values.  The  temperatures  in  the  nozzle  edge 
gradually  decrease  until  the  end  of  in-bore  time.  All  three  ceramic  nozzle  temperature  profiles 
are  higher  than  that  of  the  steel  nozzle.  The  chamber  temperature  profiles  are  unaffected  by  the 
nozzle  materials  and  are  the  same  for  each  calculation. 
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Figure  10.  Thermal  progressions  in  ceramic  nozzles  compared  with  the  steel  baseline  in  area  1. 

Figures  1 1  and  12  show  the  temperature  profiles  in  the  nozzle  throat  and  top  areas  (denoted  2 
and  3).  The  temperature  varies  from  ambient  to  their  highest,  a  stable  plateau  near  1 150  °C  for 
the  three  ceramics  which  exhibit  very  similar  temperature  profiles  while  the  steel  nozzle 
temperatures  are  about  200  °C  lower  than  those  of  the  ceramics. 

4.4  In-Depth  Temperature  Gradient 

The  large-temperature  gradient  in  the  nozzle  wall  (radial  direction)  could  cause  significant 
thermal  stresses.  In  order  to  capture  such  a  large  temperature  gradient,  a  very  fine  mesh  is 
required.  An  axisymmetric  ABAUQS  FEA  model  was  built,  see  figure  13.  In  this  model,  the 
four-node  axisymmetric  element  mesh  size  is  used  with  a  nominal  0.01 -mm  dimension.  The 
total  number  of  elements  used  in  the  model  is  about  30,000.  The  input  material  properties  used 
in  the  axisymmetric  model  are  the  same  as  those  discussed  in  section  3. 
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Figure  11.  Thermal  profiles  of  the  ceramic  nozzles  compared  with  the  steel  baseline  in  the  throat  area. 
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Figure  12.  Thermal  profiles  of  the  ceramic  nozzles  compared  with  the  steel  baseline  in  area  3. 
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Figure  13.  Axisymmetric  FEA  model  for  the  nozzle  test  fixture. 


Figure  14  shows  the  temperature  contours  for  the  steel  and  ceramic  nozzles  at  17  jus.  It  is  seen 
that  immediately  in  from  the  nozzle  surfaces,  the  temperatures  are  significantly  reduced  in  all 
cases.  At  a  distance  0.2  mm  under  the  steel  nozzle  surface  the  peak  temperature  is  775  °C,  as 
shown  in  the  figure  15.  For  the  steel  nozzle  case,  a  1030  °C  temperature  difference  is  noted  in 
0.1  mm  from  the  surface.  Similarly,  the  temperature  gradients  of  ceramic  nozzles  for  SN47, 
STK4,  and  ZR02  are  plotted  in  figure  16.  It  is  expected  that  such  large  temperature  gradients  in 
depth  will  cause  significant  thermal  stresses. 


Figure  14.  Axisymmetric  temperature  contour  plots  for  the  steel  and  ceramic  nozzles. 
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Figure  15.  Temperature  gradient  from  the  surface  in  depth  for  the  steel  nozzle. 


Figure  16.  Radial  temperature  profdes  under  the  surface  for  the  steel  and  ceramic  nozzles. 
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5.  Summary 


Three-dimensional  FEA  thermal  analysis  was  conducted  for  the  37-inm  erosion  test  fixture  for  a 
steel_4340  gun  chamber  with  a  steel  and  three  ceramic  nozzles.  Previously  computed  data  for 
the  space-dependent,  gas-temperature  and  convective-heat  transfer  coefficient  was  used  in  the 
analysis.  Temperature-dependent  heat  conductivities  and  heat  capacities  were  used  for  both  the 
gun  chamber  and  nozzles.  The  analysis  shows  that  significant  radial  temperature  gradients 
occurred  in  the  nozzles.  The  nozzle  entrance  exhibited  the  highest  surface  temperatures.  The 
three  ceramic  nozzle  temperature  profiles  show  similar  trends  and  are  approximately  200  °C 
higher  than  that  of  the  steel  nozzle,  due  to  their  relatively  low  thermal  conductivity.  The 
axisymmetric  FEA  thennal  model  was  also  built  for  in-wall  thennal  gradient  analysis.  Thennal 
profiles  under  the  nozzle  surfaces  for  all  the  materials  were  obtained  from  this  analysis  and  are  a 
necessary  input  for  the  future  thermal  stress  analysis. 
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